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Abstract 

It is shown that gravitating magnetic fields affect the evolution of curvature perturbations in 
a way that is reminiscent of a pristine non-adiabatic pressure fluctuation. The gauge-invariant 
evolution of curvature perturbations is used to constrain the magnetic power spectrum. Depend- 
ing on the essential features of the thermodynamic history of the Universe, the explicit derivation 
of the bound is modified. The theoretical uncertainty in the constraints on the magnetic energy 
spectrum is assessed by comparing the results obtained in the case of the conventional thermal 
history with the estimates stemming from less conventional (but phenomenologically allowed) 
post-inflationary evolutions. 
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1 Introduction 



It is tempting to speculate that large-scale magnetic fields are generated in the early Universe [H l2l [3] 
during a phase where the Weyl invariance of their evolution equations is broken either spontaneously 
of explicitly (see, for instance, and references therein). As a result the obtained magnetic fields will 
be necessarily tangled over typical wavelengths larger than the Hubble radius [SI EJ [TJ El [9] . In this 
situation, the presence of magnetic fields may affect the evolution of curvature inhomogeneities in a 
way that is reminiscent of what happens in the presence of a non-adiabatic pressure fluctuations. 

Suppose, indeed, that large-scale magnetic flelds are generated thanks to the ampliflcation of the 
quantum fluctuations of an Abelian gauge fleld which can be identifled with the hypercharge field 
whose modes will not be screened by thermal effects. In spite of the specific model of parametric 
amplification, the result of the process will be a stochastically distributed field whose two point function 
can be written as 

{B,{k,T)B^{p,T)) = '^pi{k)P^{k)6^'\k + p), (1.1) 

where PB{k) is the magnetic power spectrum (which may change depending upon the specific way the 
Weyl invariance is broken) and Pi{k) = [Sj — kik^/k'^) is the traceless projector. In Eq. (11.11) the 
momenta appearing in the correlator are the comoving wave-numbers. The homogeneous component 
of this configuration vanishes and, as a result, the magnetic fields will not break spatial isotropy I of 
the background geometry whose line element can be written, in the spatially fiat case and in terms of 
the conformal time coordinate r, as 

ds"" = a\T)[dT^ - dx% (1.2) 

The typical wavelengths that set the initial conditions of the CMB anisotropics § are still larger 
than the Hubble radius (i.e. kr < 1 ) at recombination, taking place for an approximate redshift 
Zrec ~ 1050 when the visibility function is peaked. When kr < 1 the evolution equations of the 
perturbations of the spatial curvature will receive contributions from three independent sources. The 
adiabatic mode typically provides the first contribution arising from the quantum fiuctuations of a 
(single) infiaton field. The second contribution may come from a collection of non-adiabatic modes 
that may be present either because the infiationary evolution is driven by more than one scalar degree 
of freedom, or because the primordial plasma contains a number of spectator fields which do not drive 
the evolution of the background but whose inhomogeneities contribute to the overall fiuctuations of 
the spatial curvature. Finally, the third contribution to the evolution of curvature fiuctuations can be 
attributed to gravitating magnetic fields. Up to now various studies addressed the interplay between 
large-scale magnetic fields and the scalar [121 [13 [HI (see also [151 [IS])) vector [TTl [TH [19] and tensor 
|19j modes of the geometry (see [2D] and [2T] for two recent reviews). As far as the scalar modes are 
concerned, it was assumed that large-scale magnetic fields are present prior to recombination (but 
after equality) and the resulting corrections to the Sachs- Wolfe plateau have been computed [T2l [13]. 
The analysis has been also recently extended to cover the Doppler region, i.e. the first, second and 
third peaks of the temperature autocorrelations fT^. The main theme of the present paper will be 

^If a background magnetic field exist, it must necessarily be oriented along a specific direction [lO]. The effective 
spatial isotropy is broken and the angular power spectrum depends on the direction of the magnetic field. As a 
consequence a number of constraints, stemming directly from the breaking of spatial isotropy, can be derived llj. 

^By this expression we mean the wavelengths that set the initial conditions of the lowest multipoles of the Boltzmann 
hierarchy. 
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to discuss possible (further) numerical constraints arising from the coupled evolution of large-scale 
magnetic fields and curvature perturbations. This analysis, when correctly performed, allows to assess 
the theoretical uncertainty of the computational scheme. In fact, the results obtained in the case of 
the conventional thermal history can be compared with the estimates stemming from less conventional 
(but phenomenologically allowed) post-inflationary evolutions. 

In the most simplistic scenario a conventional inflationary phase is followed, after a sudden 
reheating, by a radiation-dominated phase which is replaced, at equality by the standard matter- 
dominated phase. This evolutionary history allows to compute the curvature perturbations that are 
induced, for instance, by a single primordial adiabatic mode. It is also possible, in the same framework, 
to include the peculiar effect of large-scale magnetic fields whose inhomogeneities may indeed affect 
the curvature perturbations [121 1131 IIHl [K] . The result of this calculation allows the estimate of the 
Sachs- Wolfe plateau and, ultimately, the estimate of the Doppler oscillations [H]. 

If the evolutionary history of the background geometry changes in its early phases, also the effect 
of large-scale magnetic fields on the induced curvature perturbations will be different. The evolution 
of magnetized curvature perturbations is sensitive, both, to the total barotropic index and to the total 
sound speed, i.e. 



where the prime denotes a derivation with respect to r and the subscript t reminds that the barotropic 
index and the sound speed are the ones pertaining to the total fluid. 

If the rate of expansion changes in the early phases of the thermal history of the Universe (typically 
prior to big bang nucleosynthesis) also Wt and will be different and the curvature perturbations will 
have a slightly different evolution. This qualitative argument will now be corroborated by a detailed 
calculation and this is one of the themes of the present investigation. The problem then becomes 
to estimate the total curvature perturbation that will be sensitive both to the adiabatic and to the 
magnetized contributions. 

The strategy adopted in the present investigation will then be to scrutinize different thermodynamic 
histories of the Universe and to compute, in each of these cases, the resulting curvature perturbations 
for wavelengths that are still larger than the Hubble radius after matter-radiation equality. It will be 
shown that possible variations in the evolution of the total barotropic index will modify the interplay 
between the adiabatic and the magnetized components of curvature perturbations. This procedure will 
also allow to determine the theoretical error caused by selecting the most simplistic thermal history 
in comparison with its non-minimal counterparts. 

The plan of the present analysis will then be the following. In section 2 the essential theoretical 

^Notice that the most simphstic scenario is also the one contemplated by the minimal paradigm compatible with 
the three data sets that are used to infer the values of the cosmological parameters, i.e. the ACDM paradigm. These 
three data include, in general terms, the large-scale structure observations, by the CMB observations and by the type 
la supernovae observations. Before plunging into the discussion, it is appropriate to comment on the choice of the 
cosmological parameters that will be employed throughout this section. The WMAP data [22l [23l [24l [25l [26] have 
been combined, so far, with various sets of data. These data sets include the 2dF Galaxy Redshift Survey [27], the 
combination of Boomerang and ACB AR data [28l [29] , the combination of CBI and VSA data [30l [31] . Furthermore 
the WMAP 3-year data have been also combined with the Hubble Space Telescope Key Project (HSTKP) data [35] as 
well as with the Sloan Digital Sky Survey (SDSS) [33l [34] data. Finally, the WMAP 3-year data can be also usefully 
combined with the weak lensing data [351 136] and with the observations of type la supernovae (in particular the data of 
the Supernova Legacy Survey (SNLS) Jf] and the so-called Supernova "Gold Sample" (SNGS) [Ml [33). 




Pt 




(1.3) 
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tools will be introduced. A consistent gauge-invariant description will allow to follow, in one shot, 
the evolution of the curvature fluctuations and of the density contrasts on comoving orthogonal hy- 
persurfaces. In section 3 different models for the evolution of the barotropic index will be scrutinized 
and motivated. The resulting (magnetized) curvature perturbations will be computed. In section 4 
a set of bounds on the magnetic field intensity will be derived and compared. Section 5 contains the 
concluding remarks. 



2 Evolution equations 

The notations ubiquitously employed in the present script imply that the Friedmann-Lemaitre equa- 
tions in the spatially fiat case (with line element given in Eq. (11.21) ) are given by: 

= ^-^a^p,, (2.1) 
n'-n' = AiTGa\p,+p,), (2.2) 
p[ + 3n{pt+pt) = 0, (2.3) 

where Ti. = a' /a and the prime denotes, as in Eq. (11.31) a derivation with respect to r. Notice that 
H = aH where H = a/a is the Hubble expansion rate and the dot denotes a derivation with respect 
to the cosmic time coordinate t (recall that the connection between t and r is given by the differential 
relation dt = a{T)dT). By virtue of the connection between Ti. and H, a given wavelength is larger 
than the Hubble radius provided the corresponding comoving wave-number k satisfies the condition 
k/{aH) c^ikr < 1. Equations ([L3]) and (ES) imply 

3n[Wt + 1) 3 dlna 

so that = Wt iff the (total) barotropic index is constant in time. For the purposes of the present 
investigation it is practical to start with the evolution equations of the fluctuations of the geometry 
expressed in fully gauge-invariant terms. The Hamiltonian and momentum constraints can be written 

as 

V^^ - 37^(7^$ + ^') = AnGa'^idpt + 6pb), (2.5) 
V\n^ + ^') = -AnGa^ (1 + Wt) Pt^t, = Wt, (2.6) 

where Spt is the total (and gauge-invariant) density contrast of the fluid sources and 5pb{t,x) = 
52(x)/(87ra^). The fluctuations $ and ^ are the gauge-invariant Bardeen potentials [lOl HI] that 
coincide with the longitudinal fluctuations of the metric in the conformally Newtonian gauge. In 
Eq. (12. 6p dt is the three- divergence of the total peculiar velocity. At the right-hand side of Eq. 
(12.61) the three-divergence of the Poyniting vector leads to a term going as V ■ [i? x B]/a'^, where 
E ~ J /a ~ (V X B)/a is the Ohmic electric field. The conductivity a is large during most of 
the thermodynamic history of the Universe, since it is proportional to the temperature when the 
temperature is much larger than the mass of the correponding species. The contribution of the Ohmic 
Poynting vector will therefore be neglected in comparison with the other components of the magnetic 
energy-momentum tensor. Introducing the density contrast on comoving orthogonal hypersurfaces, 
i.e. 

6^(r, x) = ^-flllfR - m{l + w,)V,, (2.7) 
Pt 
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the Hamiltonian constraint (12.51) takes the pecuhar form 



V'^ = 47rG'a>tem, (2.8) 

which is (just formally) analog to the Poisson equation typical of the non-relativistic treatment of 
gravitational inhomogeneities. Another relevant pair of gauge-invariant quantities is represented by 

C = -^- ^^^P^]^P^\ (2.9) 

n = -^- ^^^^ + ^') (210) 

The variable ( [HI Il2] defined in Eq. (12. 9p . if evaluated in the gauge where the spatial curvature is 
uniform, is proportional to the total density contrast (as it can be directly checked by using Eq. (12. 3p ). 
The variable TZ [l2l HS], defined in Eq. (12.101) . if evaluated on comoving orthogonal hypersurfaces, 
coincides with the fluctuations of the spatial curvature. It is clear that the three gauge-invariant 
variables defined, respectively, in Eqs. (12. 7p . (12. 9p and (I2.10p are all related by the Hamiltonian 
constraint. Inserting Eq. (12.30 into Eq. (12. 9p and taking into account Eqs. (12. Sp . (I2.10p and (12.50 the 
Hamiltonian constraint can be expressed in one of the following two equivalent forms: 

C = ^+ i9 r 2 V (2-11) 

12nGa^Pt{l + Wt) 

^ = '^+ ',n'l V (2.12) 
3(1 + Wt) 

where Eq. (I2.12p follows from Eq. (12. lip making use of the generalized Poisson equation (12. 8p . It is 
relevant to remark that when the wavelengths are larger than the Hubble radius, at a given epoch. 

In spite of the physical differences between (, TZ and the explicit solution of the whole system in 
terms of one of these variables allows to compute the others. From the fluctuation of the covariant 
conservation equation of the total plasma the following equation can be easily obtained 

5p[ + m{5p, + 5pt) - 3(pt + PtW + (Pt + Pt)Ot = 0. (2.14) 

Using now Eq. (12.90 inside Eq. (12.140 a first-order differential equation for ( emerges and it is given 
by: 

C = 77— ^^^^nad + % B - ^. (2.15) 

pt(l + wt) 3pt(l + wt) 3 

The term 5pnad in Eq- (I2.15P accounts for the non-adiabatic pressure inhomogeneities; the following 
decomposition of the total pressure fluctuation 



r / "1^^ ' \ r / — ^st"ri. ' '^J^iiaa; (2.16) 

has been tacitly assumed in Eqs. (12.140 and (I2.15p . While the first (adiabatic) contribution gives 
the pressure fluctuation produced by the inhomogeneous density fluctuation when the specific entropy 
is constant (i.e. 5<j = 0), the second (non-adiabatic) contribution arises even if the density is un- 
perturbed (i.e. 6pt = 0) but the plasma possesses many different components (for instance, in the 
pre-recombinantion plasma, photons, baryons, neutrinos and dark matter). 
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The evolution of TZ can be obtained, in similar terms, from the equations of $ and \E', i.e. 

^" + n{<t>' + 2^') + {2n' + Ti?)^ + -V^($ - ^) = 47rGa^(5pt + ^pb), (2.17) 

3 

{d,& - \5lV'y^ - VI/) = 87rGa2(nK n^), (2.18) 
where 11^ is the anisotropic stress of the fluid and 

TT^ = (BiB' - ^bA (2. 19) 



47ra4V 3 

is the anisotropic stress of the magnetic field. Inserting Eq. ()2.10p into Eq. (12.1 7p and recalling Eq. 
fl2.16p the evolution equation of TZ can be written as 

Pt(l + wt) 3pt(l + wt) 12TcGa^{pt + pt) 47rGa^(j9t + pt) 

By now subtracting Eq. fl2.20p from Eq. (12.150 and by using Eqs. (12.120 and (12.130 . the appropriate 
equation of can be obtained after simple algebra: 

e'^ - mw^e^ = -(1 + w,)e, - m{l + wt)llt (2.21) 

where the notation 

d^dm) + d,d% = {pt + Pt)v2nt, (2.22) 

has been used. If there are coUisionless particles in the plasma (like neutrinos, after weak interactions 
have fallen out of thermal equilibrium) we will have that 

[pt + Pt) V'Ht = {p, + p.) V'n, + {p^ + p^) V'Hb. (2.23) 

Concerning Eq. (I2.23P two comments are in order. We referred the magnetic anisotropic stress to 
the photon background. This is natural since, for typical length-scales larger than the magnetic 
diffusivity scale, both quantities redshift at the same rate with the expansion of the Universe. The 
second comment involves the relevance of the anisotropic stress. As discussed in relation with the 
estimate of the CMB autocorrelation induced by magnetic fields, the magnetic anisotropic stress is of 
upmost importance at intermediate scales. In particular the interplay with the neutrino anisotropic 
stress leads to the correct initial conditions for the magnetized CMB anisotropies. However, in the 
present paper we will be mostly concerned with the largest scales and we will try to assess what could 
be the influence, on those scales, of slightly different thermal histories of the Universe. For typical 
wavelengths larger than the Hubble radius at recombination the anisotropic stress can be neglected in 
the first approximation. If needed, however, it can be included with the standard iterative procedure 
(see also HBl) where, to lowest order, the solution is the one where $ ~ v]/. 



3 Different thermal histories 

By looking at Eq. (I2.15P it is clear that, to leading order in fc^r^, the evolution of ( (and of TZ) is 
determined by three independent pieces of information: 

• the presence (or absence) of non-adiabatic pressure fluctuations; 
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• presence (or absence) of large-scale magnetic fields; 

• the specific time dependence of Wt and Cg^ 

The simplest possible situation compatible with the presence of super-Hubble magnetic fields arises 
when the Universe becomes suddenly dominated by radiation at the end of inflation. This is, inciden- 
tally, also the underlying assumption in the standard ACDM scenario. The Universe will then become 
dominated by dusty matter at a redshiftlf] 2;eq ~ 3200. If inflation was driven by a single infiaton and 
if other spectator fields with scale-invariant spectrum were absent, then it is also rather plausible to 
enforce the condition (5pnad = 0. In this situation, the exact solution of Eqs. (12.11) . (12. 2p and (12.31) 
imphes that 

1 4 

Wtioi) = —, r, C^t(«) = —, TTi (3.1) 

^ ^ 3(a + l)' ' 3(3a + 4)' ^ ' 



where 



a(r 


) = — = X^ 


+ 2x, 


X = 


T ao 








Ti ' Oeq 


ao 


^O^RO 


''"eq 


iV2- 


l)req. 



0"R0 

(3.2) 

Using Eqs. (13.11) and (13. 2p inside Eq. (I2.15p and changing the variable from r to a the following 
expression can be obtained: 

= — 71x2^ ^7 = — TT ^b{t,x) = -— . (3.3) 

da (3a + 4)2 p^{t) p^{t) 

In Eq. (13. 3p the magnetic energy density has been expressed in units of the photon energy density and 
denotes the photon fraction in the post-inflationary radiation background. After neutrinos decouple 
for temperatures of the order of the MeV the photon fraction will be related with the neutrino fraction 
as Rj = 1 — Ri, where R^, = 0.405 (for three neutrino families). Direct integration of Eq. (13. 3p gives 
the sought result, namely 

(3,4) 

The constant C*(^) stands for the adiabatic mode produced (for instance during a phase of conventional 
inflation). After matter-radiation equality (but before recombination) to leading order in fc^r^ we will 
have, from Eq. (12.110 . that 

Rf^sik) f 3aiec \ , ^^/j 2 2 \ ^cq + 1 _ o nv ^O^MO 



4 \3arec + 4/ Zj-ec +1 V 0.134 / 

This result can be easily interpreted from the physical point of view. When the rate of expansion 
increases from radiation to matter the barotropic index (and the sound speed) are both decreasing 
from 1/3 to 0. Thus the overall effect on the source term for the evolution of curvature perturbations 
implies that the magnetic contribution tends to cancel the contribution of the adiabatic mode. 



■''When not otherwise stated it will be assumed that HqUmq — 0.1326, where Hmo denotes the present fraction in 
dusty matter. This value for the total critical fraction of matter emerges, in the context of the ACDM paradigm, when 
the WMAP data [22l [23] are combined with all the remaining cosmological data sets of different origin (excluding weak 
lensing data). 
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In the context of this reference scenario when the Universe only passes, after inflation, from radia- 
tion to matter also and can be computed, respectively, from Eqs. (12. 8p and fl2.6p . It suffices to 
extract ^ from C, and then use the aforementioned constraints. Going to Fourier space result of this 
manipulation will be 



9t{k,T) 



V2a\fa 

k^Ti 



C(A:)>Vi(«)--i?^(^B(A;)W2(«) 



a 



2a2 



Uk) + ^R'yflBik) 



2a' 



(3a + 4)y/a + 1 



(3.6) 
, (3.7) 



where 



Wi{a) 



15y/a + 1 
2 



{16[v^ - 1] + a[a{9a + 2) - 8]}, 



5^/a + l 



{16[1 - V^^TT] + a[8 + a{a - 2)]}. 



(3.8) 



Thanks to the success of big-bang nucleosynthesis, it is rather plausible to imagine that the Universe 
was already dominated by radiation for temperatures as large and few MeV. However, prior to that 
time there are no direct probes of the expansion rate of the Universe. In a somehow indirect way this 
could be achieved through the study of the cosmic background of gravitational waves. At the moment, 
for instance, there is no compelling evidence of the fact that, after inflation, the Universe became 
suddenly dominated by radiation. It is, on the contrary, rather reasonable that, before settling on a 
radiation-dominated stage of expansion the Universe passed through a phase where, for instance, the 
rate of expansion was smaller than the Hubble rate during radiation. This occurrence is realized when 
the total barotropic index of the sources driving the geometry satisfles 1/3 < z/^t < 1- In this case the 
sources are said to be stiffer than radiation [47J. The stiffest equation of state we can imagine is the one 
where the barotropic index is 1 which implies that the sound speed coincides with the speed of light. 
The case of speed of sound equal to speed of light is the one contemplated by the so-called Zeldovich 
model (see [50] and references therein). A similar kind of evolution arises in quintessential inflationary 
models where the inflaton and the quintessence field are unified in a single degree of freedom. In this 
case (as in the more general case of stiffer phases [17]) the spectrum of relic gravitons presents a sharp 
increase [51] and a peak [52] which are potentially accessible to direct observations. In the context of 
quintessential infiation the stiff epoch is effectively dominated by the kinetic energy of the inflaton. 
Stiff phases also arise in brane-world scenarios [S3] where the calculation of the graviton spectrum 
mirrors the four- dimensional case [STl [521. 



The transition from a stiff epoch to a radiation-dominated phase of expansion can be understood 
in terms of the following exact solution: 



a 



Or 



r 

2/ = — , 



where 




Mp 



A 



'Mp 



(3.9) 



(3.10) 



Equations (13.91) and (I3.10p are solutions of Eqs. (12. ip . (12. 2p and (12. 3p when the sources of the geometry 
are given in terms of a mixture of radiation and of a stiff fluid with Ps = Ps- Equations (13.90 and (13.100 
parametrize the physical situation where, at r^, some (small) amount of radiation is present in compar- 
ison with the stiff contribution (i.e. pr* ^ ps*)- In the context of quintessential inflationary models 
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the initial amount of radiation comes indeed from quantum fluctuations [55] and, consequently, 
Pr* ~ H^. This occurence implies, in turn, that the hierarchy between the initial scale factor at r* 
and the scale factor at the onset of radiation (i.e. is determined by M-p/H^). Since if* < 10~^Mp, 
the duration of the stiff phase is correctly bounded because Hy- > H\,hn- According to Eq. (12. 4p . the 
barotropic index and the sound speed are: 



1^2 + 3 



2a^ + 9 



3a2 + l' 3(2a2 + 3)- 

Using therefore Eq. (13. lip inside Eq. (12.150 and changing variable from r to a: 

da ~ (2a2 + 3)2' 
Integrating Eq. (I3.12p between a = and ttf we do get 



(3.ii: 



(3.12) 



ak,T) = uk) + 



(3.13) 



2{2al + 3) 

where af ^ 1 and it is defined deep in the radiation epoch. This result implies that the magnetic 
contribution enhances (rather than canceling) the adiabatic term C*(^) of inflationary origin. 

As in the case of the radiation-matter transition the total density contrast and the total peculiar 
velocity can be computed. In Fourier space the result is: 



6v^ 



(3.14) 



9,ik,T) = -k\ 
R^QB{k) 



a 



2WM 



{2a^ + 3)V(y^ + 1 



r>V4(al 



where 



1] 



(3.15) 



2a' 



2{2a' + 3) 



(3.16) 



This result can be generalized when the stiff phase is parametrized in terms of a source = 
with 1 < (5 < 1/3. Repeating the same steps outlined above we will have that the dependence of the 
barotropic index and of the sound speed upon the scale factor are given, in this case, by: 



1 /3<5 + a 



35-1 



3 V 1 + «35-i 



35-1 



9(5(5+ 1) + 4a 
3[3(5 + l)+ 403-5-1] 



(3.17) 



The equation for and its related solution can the be written as 

da " ^ ^ [3(5 + l)+4a3<5-i]2 ' 



C(fc,r)=C(A;) + i?^fiB(fc)- 



a 



35-1 



4af -1 + 3(5 + 1) 



(3.18) 
(3.19) 
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As in the other two cases em and 6^ can be easily computed: 

^2 2 a'{S-l)/2 



Ot{k,T) = 



952 - 


1 VI + a^^~^ 


k 


^(35+l)/2 



V36 



a 



3(5-1 



3(5-1 



+ 1 



(3.20) 



Va=^^-^ + lW5(a,5) 



2(35-1) a(35+5)/2 



35-1 



+ 1 W6(a,5) 



4a35-i + 3(5 + 1) a(35+5)/2 2(35-1) 



(3.21) 



where 



K") 4s + 3(5 + 1) 7-3* , 

^ ^.S 2(35-1) 



5+3* 
s{a) 5 2{3*-l) 



5+1)3/2 

ds. 



[S + 1)3/2 



3(5-1 



(3.22) 



For each value of 6 the 



The upper limit of integration in Eqs. fl3.22p corresponds to s{a) 
above integrals can be evaluated in explicit terms. 

Suppose now to investigate a three stage model where the Universe expanded, initially, at a rate 
that was slower than radiation then got trapped in a radiation phase which turned subsequently into 
an epoch dominated by dusty matter. In this case we will have that the barotropic index and the 
sound speed are given, respectively, by: 



Wt 



3 + 



31 + a^ + ba^' 
18 + 4^2 



Pm* I Ps*_ 
Ps* \PR* 



3/2 



18 + 12a2 + 9^^3- 

In this case the evolution equation of curvature perturbations is given by 

da ^ ^(3a36 + 4a2 + 6)2■ 



(3.23) 
(3.24) 



(3.25) 



The (approximate) solution of Eq. (13.251) in the matter-dominated epoch reads: 

R^riBik) 



C{k,r)^Uk) + 



3arec + 4' 



(3.26) 



The result of Eq. (13.261) shows an extra suppression of a factor a^-ec — 1/3 in comparison with 
Eq. (13. 5p . This suppression can be easily understood. The evolution described by Eqs. (13.231) and 
( ]3.24p implies that, right after inflation, the Universe expands slower than radiation. When radiation 
kicks in, the curvature perturbations are enhanced. Later on, when dusty matter becomes dominant 
(around matter-radiation equality) the magnetized contribution tends, again, to cancel the pre-existing 
adiabatic mode. The net result of the initial increase and of the subsequent decrease of ( is given in 
Eq. (I3.26P and it is a bit different from the result of Eq. (13. 5p where the intermediate stiff phase was 
absent. The following conclusions can then be drawn: 

• if the Universe is dominated by radiation from the end of inflation (as it happens in the case of 
the standard evolution) the magnetized contribution and the adiabatic mode have opposite sign; 
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• if, prior to the electroweak epoch, the thermodynamic history of the Universe deviates from a 
radiation epoch the magnetized contribution is more suppressed in comparison with the standard 
case. 

In the framework of conventional inflationary models this type of deviation from the standard sce- 
nario is pretty general. By general we mean that, within the current bounds on the expansion of 
the Universe arising from big-bang nucleosynthesis we can just imagine drastic deviation from the 
radiation-dominated evolution between the end of inflation and, say, neutrino decoupling. In fact, 
deviation from a radiation-dominated evolution after (or right before) big bang nucleosyntheis will 
necessarily jeopardize the production of the abundances of the light elements. 

Up to now it has been assumed that the fluids of the primeval plasma do not exchange energy. Let 
us now address this interesting situation that could arise, for instance, in the course of reheating. Let 
us therefore imagine that a matter fluid (which can model the coherent oscillations of the inflaton) 
decays into massless particles. This dynamics can be parametrized by the following system |^: 

ij = -47rG(^Pr + Pm), (3.27) 
= -SHp^, - Lpm, (3.28) 
Pr = -4/fp, + Tp^. (3.29) 

The Universe, initially dominated (right after inflation) by dusty matter, becomes rapidly dominated by 
radiation at a rate that is controlled by F. It should be immediately noticed that, by summing up Eqs. 
f l3.28p and fl3.29p the total energy density pt = Pr+Pm is covariantly conserved, i.e. pt + 3H{pt+pt) = 0. 
Equations fl3.28p and (13.291) can be approximately solved: 

Pm(t) = Pr(ti)(^) V^(*"*'), pr(t) ^ PriU)[^)\ (3.30) 

where — T \ As a consequence of the relations of Eq. fl3.30p we can also say that 

Pr H^\aJ \adJ 

In this situation the total barotropic index and the total sound speed are slightly modified thanks to 
the physical differences of the system. In particular we will have that 

c2 = 1 Pi E PllL (-3 32) 

34pr + 3p,„ 3H4p, + 3p^' ^ ' ' 

From Eq. (l3:32|) is clear that 4 ~^ 1/3 when F > /f. 

The evolution equations of curvature perturbations to be solved read: 



C = ^— 5pnad + ^ Ult - ^"I^Pb, (3.33) 

Pt + Pt Pt + PtV 3/ 

Cm = -Cm-r(7^C, (3.34) 



^Here we passed from the conformal to the cosmic time parametrization. See the comments after Eqs. p.ip . (|2.2p 
and (ED). 
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where g = H/{T + 3iJ). Let us now concentrate, as discussed above, on the conventional situation 
where the non-adiabatic modes are totally absent. In this case, initially, (Jpnad (Cr — Cm) — and 
Cr = Cm — C*- Using Eq. (13.320 inside Eq. ( 13.33^ we will simply have that 

/ _ H6pBpm .oocrN 

Recalhng that, after td, T ^ H and pm is exponentially suppressed, Eq. (I3.35P can be written as 

^ = -MR^e-^, x = Tt. (3.36) 
ax lb 

Integrating once Eq. (13.361) the result will be 

ak,y) = Uk) - ^^^^ r v^e-^ dx. (3.37) 
16 Ji 

The integral mentioned above can be done introducing the error function, i.e. 

Taking the limit y oo, we then have 

C(M)^C.W-^(^i^^^^^^^). (3.39) 



that is, numerically. 



ak,t) ~ Uk) - 0.079 R,QM ~ Uk) - ^h^M^, (3.40) 



4 Bounds on different thermal histories 

The two-point correlation function of curvature perturbations will receive contributions, in general, 
from the adiabatic mode, from one (or more) non-adiabatic modes, and from the magnetic field. As 
argued in the previous sections, the simplest situation is the one where only the adiabatic mode is 
present together with the magnetized contribution that we ought to constrain. This choice is motivated 
by the experimental evidence that the Doppler peak structure of the temperature autocorrelations 
strongly suggests that, after equality, the large-scale curvature perturbations were, predominantly 
adiabatic [221 [23l [211 [251 [26] . 

Consider therefore, the curvature perturbations present at recombination in the framework of the 
different thermal histories discussed in the previous section. The curvature perturbations can be 
written, in Fourier space, as0 

Ch{k, ttrec) = C*(^) + Th{arcc, /io^Mo)fiB (A^) , (4.1) 

where the subscript h refers to each different thermal history and where C*(^) represents the adiabatic 
contribution normalized to the large-scale value of the (ordinary) Sachs- Wolfe contribution to the 

^The quantity Th(arec, hQiluo) is the generalized transfer function that may change depending upon the specific 
thermal history labeled by the subscript h. 
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temperature autocorrelations. The total (scalar) power spectrum Vh{k), i.e. the Fourier transform of 
the two-point function computed at different spatial locations but at the same time is defined as 

/SI Tl nT" 
d\nkVh{k)^—, r = \x-y\. (4.2) 

From Eq. (14. ip and within the conventions summarized by Eq. (14. 2 p we then have that 

(a(fc,T)a(p,r) = ^p,(A;)5(3)(fc+p-) (4.3) 

where 

Vh{k) = Vcik) + T2(«,,„ hlnMo)Vn{k) + 2Th{a,,,, hlVlMo)^V^{k)^Vn{k) cos 7, (4.4) 

where 7 is the correlation angle between the adiabatic mode and the magnetized contribution. In Eq. 
(14.41) Vc^ is the power spectrum of the adiabatic contribution defined as 

/ k \ "c-i 2 X 1 

n = A r , A = — ^A(fcp) = 2.95 X 10-9A(A;p). (4.5) 
\ Kp / yi ^^^^ 

In Eq. (14. 5 p is the spectral index of the adiabatic mode and Ac_ is the amplitude referred to 
the pivot scale kp. Following the WMAP conventions the value of the pivot scale will be fixed as 
fcp = 0.002 Mpc~^. The numerical factor appearing in the expression of Ac^ has been obtained by 
taking consistently Tcmb = 2.725 x 10^ (expressed in units of ^K). In the absence of any other 
contributions the WMAP 3-year data implj]^ ~ 0.947 and A{kp) ^ 0.815 when combined with the 
remaining cosmological data sets, i.e. supernovae and large scale structur^. In Eq. (14. 5 p the magnetic 
part of the correlation function is expressed as |T3l [H] 



,2 / \ , , 4(6 -e)(27r)2^ 



= , m = (4.6) 

where 2e the spectral index of the magnetic energy density fluctuations and A;l is the magnetic pivot 
scale that will be defined in a moment. In Eq. (14.60 we have that 

nBL = ^ = r.56xlO-^f^y. (4.7) 



Inp^ VnG 

In Eq. (14. 7p i?L is the value of the magnetic field smoothed, through a Gaussian window function, 
over a typical comoving length L which is related to the magnetic pivot scale as L = 2n /ki^. In 
what follows the fiducial value = Mpc~^ will be consistently adopted. It is relevant to appreciate 
that i?L represents the smoothed magnetic field red-shifted to the present epoch. This convention 
is a bit confusing but we will follow it since it is a common practice in the field. The source of 
the possible confusion is, in short, the following. The field -Bl is the value of the magnetic field at 
recombination redshifted to the present epoch and assuming magnetic flux conservation which is well 
justified since the value of the conductivity is always rather large. However, Bl is not the present 
value of the magnetic field intensity observed in galaxies and galaxy clusters. The rationale for this 



®In this analysis the values of the cosmological parameters will be fixed to the best fit values of the WMAP data 
combined with the other sets of cosmological observables. This choice is not crucial since the values of the cosmological 
parameters are purely illustrative. The essential features of the present analysis are unchanged if the best fit values of 
the WMAP data alone (or partially combined with the other data sets) are consistently assumed. 

^The data [32] and [351 [3S] are not included in the joined analysis. 
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statement is that during the gravitational collapse of protogalaxies the magnetic field intensity 5l will 
be amplified by compressional amplification and, probably, also by some dynamo action. Assuming 
just compressional amplification (which is the most certain aspect of the dynamics) the amplification 
factor may be of the order 10^ or even 10^. 

Notice that, from Eq. (14. 6p the nearly scale-invariant limit is achieved for e < 1. Furthermore, 
as discussed elsewhere, the power spectrum of the magnetic energy density, being quartic in the field 
intensities, is computed in terms of the appropriate convolutions that have been evaluated, to get to Eq. 
(14.61) . for e < 1. It should be borne in mind that, in the opposite case (i.e. e > 1) the magnetic energy 
density has a violet spectrum. This implies that, at very large length-scales the gravitational effect 
will be negligible and that, furthermore, the most significant bounds on the intensity of the magnetic 
field will come from comoving momenta close to the diffusion scale. Notice that the nearly scale- 
invariant spectrum is rather suggestive also because it would imply a nearly scale-invariant spectrum 
of large-scale magnetic fields from galaxies, to clusters, to super-clusters. This observation is, today, 
beyond our observational capabilities. However, by looking at recent determinations of magnetic fields 
in clusters and in some typical supercluster, it is indeed tempting to speculate that the resulting power 
spectrum of the magnetic energy density is nearly scale-invariant. 

The correlation function of curvature perturbations enters directly the determination of the or- 
dinary Sachs- Wolfe contribution which is the dominant source of temperature anisotropics at large 
angular scales. Let us require that the magnetized contribution is always much smaller than the 
adiabatic contribution. So, we will assume, in a rather conservative approach, that the magnetized 
contribution is always smaller (by a factor rj) than the adiabatic contribution. The parameter t] can 
be O{10~^), in a conservative approach. Following this strategy the magnetic field intensity can be 
bounded in each thermal history of the Universe. Consider, for instance, the plots appearing in Fig. 

where, on the vertical axis, the base 10 logarithm of the smoothed magnetic field is reported in 
units of nG. In the plot at the left the three different lines correspond to the cases of the three different 
histories discussed in the previous sections. These three different histories will imply three different 
functions Th with h = 1,2,3 and with, in general Ti 7^ T2 7^ T3. At late times the functions can be 
approximated as 



The function Ti arises in the simplest case, i.e. when the Universe, after a phase of conventional 
inflation, passes from a radiation-dominated epoch to the matter stage. In Ti the reheating is assumed 
to be a sudden process. The function T2 describes the situation when there is an intermediate stiff phase 
expanding at a rate slower than radiation. Finally the function T3 arises when there is a prolonged 
reheating and then the Universe is, subsequently, dominated by radiation and matter. In Fig. [T] (plot 
at the left) we illustrate the bounds as a function of the magnetic spectral index. The region below 
each of the curves implies that the contribution of the magnetic fleld to the two-point function is 
smaller than 10~^. 

The less restrictive case is, surprisingly enough, T2. The rationale for this occurrence can be simply 
understood in terms of the results reported in the previous section. In the transition from the stiff 




(4.10) 




(4.9) 
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Figure 1: The bounds on the smoothed magnetic field intensity are illustrated as a function of the 
magnetic spectral index (plot at the left) and as a function of the critical fraction in dusty matter 
(plot at the right). The pivot scales kp and /cl are fixed, respectively, to 0.002 Mpc~^ and Mpc~^. 

phase to the radiation phase the magnetized component gets reduced while in the subsequent transition 
from radiation to matter the magnetized contribution gets enhanced. Since the two effects tend to 
cancel, the net result will be an overall suppression of the magnetized contribution. This dynamical 
suppression allows the magnetic field amplitude to be larger than in the other two cases (i.e. Ti and 
T3) where there is no "destructive" interference between the two subsequent effects. 

In the plot at the right of Fig. [1] this pattern is confirmed with an extra piece of information: if 
the matter fraction increases the bound become looser in the case of T3 and a bit tighter in the case 
of Ti and T2. In fact, in Fig. [T] (plot at the right) the magnetic spectral index has been kept fixed 
while the matter fraction is allowed to move from the fiducial value h^Quo = 0.1326 which is the one 
assumed in the left plot of Fig. [H The increase (or decrease) of the critical fraction of dusty matter 
simply means, physically, that the recombination epoch is either delayed or anticipated. 

It should also be remarked that the case of prolonged reheating is the most constraining. This 
occurrence can be also simply understood from the considerations of the previous section. In the case 
of prolonged reheating c^^ 7^ 1/3. So during this phase the curvature perturbations decrease and this 
decrease will interfere constructively with the usual radiation matter transition since both contribution 
have the same sign. In Fig. [T]the value of the correlation angle has been fixed to 7r/2 implying that 
the adiabatic and the magnetized mode are uncorrelated. Again we may relax this assumption and 
allow 7 to vary while all the other parameters are fixed to their fiducial values. This has been done 
in Fig. [21 As expected the bounds is tighter in the case when the two components are correlated and 
looser when the two components are anti-correlated. In Fig. [2] (plot at the right) the variation of the 
adiabatic spectral index is illustrated always in the case where the two components are uncorrelated. 
It is now relevant to remark that while P({k) slightly decreases as k increases, the opposite is true 
for Pni^k). It is therefore clear that, given this situation, the most constraining wavenumbers are the 
largest, i.e. the small-scale behaviour of the spectrum is the one that leads to the most stringent 
bounds. This consideration has been taken into account in deriving the bounds expressed by Figs. [T] 
and O In particular, it has been assumed that the diffusion scale is the one roughly associated with 
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n^= 0.947, e=0.1 , Q.^^=0A326 




e=0.1 , y= 71/2, =0.1326 



0.9 0.92 0.94 0.9 



1.02 1.04 1.06 1.0 



Figure 2: The bounds on the smoothed magnetic field intensity are illustrated as a function of the 
correlation angle 7 (plot at the left) and as a function of the adiabatic spectral index (plot at the 
right). All the other parameters are kept fixed to their fiducial values. 




Figure 3: The potential differences in the estimate of the diffusion scale are illustrated. In the plot at 
the left, the bunch of curves labeled by B corresponds to the estimate of the diffusion scale given in 
Eq. (14.121) while the curves labeled by A are the same illustrated in Fig. [1] (plot at the left). In the 
plot at the right the two classes of curves have illustrate the same phenomenon but in terms of the 
dependence on and should be compared with the right plot appearing in Fig. O 
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the Silk wavenumber, i.e. 

= 9.63 X 10^3 -^-^ . 4.11 

AjDTrec V0.023/ V0.134/ V Zrec / 

There might be shghtly different choices for the diffusion scale. For instance it has been noticed in the 
past that the diffusion scale of the magnetic fields should be related with the induced velocity of Alfven 
waves. In this case the diffusion length-scale will be slightly smaller than the Silk scale by a factor 
which is essentially proportional, in our notations, to (Bi^/nG). In this case the diffusion wavenumber 
will be given by [19]: 

/ R \ -2/(e+2) 

kj,c^il.7xl0Y^(^) Mpc-^ (4.12) 



.nG. 

Figure [3] illustrates the situation described by Eq. (14.12^ . It is clear that the patterns of the dif- 
ferent thermal histories remain the same. However, the bound is improved by roughly one order of 
magnitude. 



5 Concluding remarks 

There are various lessons that can be drawn from the exercise reported in this paper. The main 
question we addressed has been the possible influence of slight variations in the thermal history of the 
Universe on the curvature perturbations induced by a magnetic field present for typical scales that are 
larger than the Hubble radius at recombination. Various inflationary mechanisms for the generation of 
cosmic magnetic field fall in this situation. It has been found that, indeed, the thermal history of the 
Universe may either suppress or increase the relative weight of magnetized curvature perturbations. 
The suppression occurs when, during the dynamical evolution, there is a sort of destructive interference. 
This phenomenon takes place, amusingly enough, when there was a phase, in the life of the Universe, 
when the rate of expansion was slower than the one of radiation. This kind of evolution must occur prior 
to the onset of big-bang nucleosynthesis and explicit examples have been provided. Also the opposite 
phenomenon can be realized, i.e. a sort of constructive interference when the overall contribution of 
the magnetized curvature perturbations gets enhanced in comparison with the adiabatic contribution 
present at the end of inflation. This situation takes place if the reheating phase is sufficiently prolonged 
in comparison with the sudden reheating approximation. 

These considerations can certainly be developed along different directions. It is interesting to 
remark that the nature of the bounds on the magnetic field intensity may change quantitatively and 
qualitatively. Different thermal histories imply that the constraints on the magnetic field intensities 
may vary, generally speaking, of one order of magnitude. In the case when the Universe contains after 
inflation (but before radiation) a stiff phase the present value of the smoothed magnetic field must be 
qualitatively!^ i?L < 1.5 x lO^^G while in the context of prolonged reheating it must be as small as 
1.6 X 10-* G. 

The presented considerations are not totally insensitive on the estimate of the diffusive wavenumber. 
In particular, it occurs that while experimental data favour a slightly red spectrum of curvature 

^°It should be stressed, as discussed in section 4, that Bj^ is the intensity of the magnetic field at recombination 
redshifted to the present epoch. This field is by no means equal to the present magnetic field. So the bounds on i?L are, 
really and truly, bounds on a primordial magnetized background. The present magnetic field of galaxies and clusters 
is certainly related to Bi^. In the simplest case, compressional amplification may turn a 0.1 nG field (coherent over a 
comoving scale of the order of the Mpc) into the /iG field observed in galaxies. 
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perturbations (i.e. an adiabatic spectral index slightly smaller than 1), the magnetized contribution is 
typically nearly scale-invariant but slightly blue. In this situation the most constraining wavenumbers 
are the ones close to diffusion scale. If the diffusivity length-scale is slightly smaller than the Silk scale, 
then the overall bounds on the smoothed magnetic field intensity range, depending on the details of 
the specific thermal history, from < nG to Sl < 0.1 nG. It can be discussed if the Alfven velocity 
should enter the thermal diffusivity scale when the magnetic field is fully inhomogeneous. Indeed 
Alfven waves are only excited when there is a background magnetic field in the game. This is not 
the case when the magnetic field does not break spatial isotropy which the case investigated in the 
present paper and which is also, in our opinion, the most realistic one (see the introduction). With 
these caveats, however, it is important to notice that different thermal histories may either relax or 
improve the bounds on the magnetic field intensity by one order of magnitude. 
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